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Abstract 

We present results for the three-loop universal anomalous dimension 7«ni(j) of Wilson twist- 
2 operators in the Af = 4 Supersymmetric Yang-Mills model. These expressions are obtained 
by extracting the most complicated contributions from the three loop non-singlet anomalous 
dimensions in QCD which were calculated recently. Their singularities at j = 1 coincide with 
the predictions obtained from the BFKL equation for = 4 SYM in the next-to-leading order. 
The asymptotics of 7«nj(j) at large j is in an agreement with the expectations based on an 
interpolation between weak and strong coupling regimes in the framework of the AdS/CFT 
correspondence. 



Introduction 

The anomalous dimensions (AD) of the twist-two Wilson operators govern the Bjorken scal- 
ing violation for parton distributions in a framework of Quantum Chromodynamics (QCD) PP. 
These quantities are expressed through the Mellin transformation 



of the splitting kernels Wb^a{x) for the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
equation ^ which relates the parton densities /a(a^,Q^) (hereafter a = X, g, (p for the spinor, 
vector and scalar particles, respectively) with different values of as follows 



The anomalous dimensions and splitting kernels in QCD are well known up to the next-to- 
leading order (NLO) of the perturbation theory p. 

The QCD expressions for AD can be transformed to the case of the Af = 1 Supersymmetric 
Yang-Mills theories (SYM) if one will use for the Casimir operators CA,CF,Tf the following 
values Ca = Cp = Nc, Tf = Nc/2 (the last substitution follows from the fact, that each 
gluino Xi being a Majorana particle gives a half of the contribution for the Dirac spinor). 
For extended supersymmetric theories the anomalous dimensions cannot be obtained in this 
simple way, because additional contributions coming from scalar particles should be also taken 
into account Recently these anomalous dimensions were calculated in the next-to-leading 
approximation for the A/" = 4 Supersymmetric Yang-Mills theory. 

It turns out, that the expressions for eigenvalues of the AD matrix in the J\f = 4 SYM can be 
derived directly from the QCD anomalous dimensions without tedious calculations by using a 
number of plausible arguments. The method elaborated in Ref. [5j for this purpose is based on 
special properties of solutions of the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation [7] in 
this model and a new relation between the BFKL and DGLAP equations (see jS]). In the NLO 
approximation this method gives the correct results for AD eigenvalues, which was checked 
by direct calculations in Its properties will be reviewed below only shortly and a more 
extended discussion can be found in 

Next-to-next-to-leading order (NNLO) corrections to AD in QCD were calculated recently [S]* 
in the nonsinglet case. Using these results and the method of Ref. [5. we derive in this pa- 
per the eigenvalues of the anomalous dimension matrix for the A/" = 4 SYM in the NNLO 
approximation. 

Evolution equation in AT = 4 SYM 

The reason to investigate the BFKL and DGLAP equations in the case of supersymmetric 
theories is based on a common belief, that the high symmetry may significantly simplify the 

*see also Ref. [S] for the singlet case 
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structure of these equations. Indeed, it was found in the leading logarithmic approximation 
(LLA) jini, that the so-called quasi-partonic operators in A/" = 1 SYM are unified in super- 
multiplets with anomalous dimensions obtained from universal anomalous dimensions 'jumij) 
by shifting its arguments by an integer number. Further, the anomalous dimension matri- 
ces for twist-2 operators are fixed by the superconformal invariance JUI- Calculations in the 
maximally extended A/" = 4 SYM, where the coupling constant is not renormalized, give even 
more remarkable results. Namely, it turns out, that here all twist-2 operators enter in the 
same multiplet, their anomalous dimension matrix is fixed completely by the super-conformal 
invariance and its universal anomalous dimension in LLA is proportional to \E'(j — 1) — ^E'(l), 
which means, that the evolution equations for the matrix elements of quasi-partonic operators 
in the multicolour limit A'c — >■ oo are equivalent to the Schrodinger equation for an integrable 
Heisenberg spin model ^] . In QCD the integrability remains only in a small sector of 
the quasi-partonic operators [12]. In the case of A/" = 4 SYM the equations for other sets of 
operators are also integrable Evolution equations for quasi-partonic operators are 

written in an explicitly super-conformal form in Ref. [T7] . 

Similar results related to the integrability of the multi-colour QCD were obtained earlier in 
the Regge limit Moreover, it was shown that in the A/" = 4 SYM there is a deep relation 
between BFKL and DGLAP evolution equations. Namely, the j-plane singularities of AD of the 
Wilson twist-2 operators in this case can be obtained from the eigenvalues of the BFKL kernel 
by their analytic continuation. The NLO calculations in A/" = 4 SYM demonstrated jS], that 
some of these relations are valid also in higher orders of perturbation theory. In particular, 
the BFKL equation has the property of the hermitian separability, the linear combinations 
of the multiplicatively renormalized operators do not depend on the coupling constant, the 
eigenvalues of the anomalous dimension matrix are expressed in terms of the universal function 
JuniU) which can be obtained also from the BFKL equation The results for 7unj(j) were 
checked by the direct calculations in Ref. ^ 

In the A/" = 4 SYM theory we have the following field content: one gluon g, four Majo- 
rana fermions A and three complex scalars 0. All particles belong to the adjoint representation 
of the gauge group SU{Nc). This model possesses an internal S'f/(4) symmetry. In the Af = 4 
SYM theory one can introduce the following colour and SU{4) singlet local Wilson twist-2 
operators [ZOl CUl 1211 122 : 









(1) 


Qg 


■■,H 




(2) 








(3) 


qX 


■■,H 




(4) 








(5) 



where T>^ are covariant derivatives. The spinors Aj and field tensor Gp^ describe gluinos and 
gluons, respectively, and (pr are the complex scalar fields appearing in the A/" = 4 supersym- 
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metric model. Indices i = 1, ■ ■ ■ ,4 and r = 1, ■ ■ ■ , 6 refer to SU (4) and 5*0(6) ~ SU (4) groups 
of inner symmetry, respectively. The symbol S implies a symmetrization of each tensor in the 
Lorentz indices /ii, and a subtraction of its traces. The anomalous dimension matrices 
can be written for unpolarized and polarized cases, respectively, as follows 



Tunpol 



799 IgX 730 
IXg 7AA 7A0 
l<Pg l<p\ 14>4> 



7poi 



733 73A 

7A3 7aa 



(6) 



Note, that in the super-multiplet of twist-2 operators there are also operators with fermion 
quantum numbers and operators anti-symmetric in two Lorentz indices [TUl I22j . For the case 
J\f = 4 the multiplicatively renormalized operators were found in an explicit way and their 
universality properties for all orders of perturbation theory were formulated in Refs. 13]. 

After their diagonalization, the new unpolarized 7 and polarized 7 AD matrices have the 
following form 



7 



V Sunpol^ 



7+(j) 7+o(j) 7+-(j) 
7o+(j) 7o(j) 7o-(j) 
7-+(j) 7-0 (j) 7-(j) 



7 



V- 



7poi 



V 



7+(j) 7+-(j) 
7-+(j) 7-(j) 



(7) 



which corresponds to AD matrices for multiplicatively renormalizable linear combinations of 
operators (H}-©. Here, the matrices V, V~^, V and were calculated in [S] and in LO we 
have 'jimU) = O5 limU) = for /,m = +,0, — . In NLO the AD matrices become triangle |lj 
due to superconformal invariance breaking [SB], similar to the case of = 1 SYM [2^. The 
eigenvalues 7z(j) and 7/(j) govern the power-like violation of the Bjorken scaling for the parton 
distributions. 

Due to the fact that all twist-2 operators belong to the same supermultiplet the anomalous 
dimensions 7i(j) and 7z(j) (/ = +,0, — ) have the properties [12] l5j 

7+(j) = 7+(j - 1) = 7o(j - 2) = 7_(j - 3) = 7-(j - 4) = 7n™(j), (8) 

where juniij) is the universal anomalous dimension. 

Method of obtaining AD eigenvalues in AT = 4 SYM 

As it was already pointed out in Introduction, the universal anomalous dimension can be 
extracted directly from the QCD results without finding the scalar particle contribution. This 
possibility is based on deep relation between DGLAP and BFKL dynamics in the A/" = 4 

SYM 

To begin with, the eigenvalues of the BFKL kernel turn out to be analytic functions of 
the conformal spin \n\ at least in two first orders of perturbation theory [H]. Further, in the 
framework of the DR-scheme [21] one can obtain from the BFKL equation (see [3]), that there 
is no mixing among the special functions of different transcendent ality levels i \ i.e. all special 

tNote that similar arguments were used also in [221 to obtain analytic results for contributions of some 
complicated massive Feynman diagrams without direct calculations. 
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functions at the NLO correction contain only sums of the terms ~ = 3). More precisely, 

if we introduce the transcendentality level for the eigenvalues of integral kernels of the BFKL 
equations as functions of 7 and appearing in the perturbation theory in an accordance with the 
complexity of the terms in the corresponding sums 

$-1/7, ~ C(2) ~ l/7^ $" ~ /3" ~ C(3) ~ l/7^ 

then for the BFKL kernel in the leading order (LO) and in NLO the corresponding levels are 
i = 1 and i = 3, respectively. 

Because in A/" = 4 SYM there is a relation between the BFKL and DGLAP equations 
(see ISlEl)! the similar properties should be valid for the anomalous dimensions themselves, 
i.e. the basic functions 'Juniij)^ luniU) luniU) assumed to be of the types ~ 1/f with 
the levels i = 1, i = 3 and i = 5, respectively. An exception could be for the terms appearing 
at a given order from previous orders of the perturbation theory. Such contributions could be 
generated and/or removed by an approximate finite renormalization of the coupling constant. 
But these terms do not appear in the DR-scheme. 

It is known, that at the LO and NLO approximations the most complicated contributions 
(with i = 1 and i = 3, respectively) are the same for all LO and NLO anomalous dimensions 
in QCD [I| and for the LO and NLO scalar-scalar anomalous dimensions |4,. This property 
allows one to find the universal anomalous dimensions 'Juniij) and 7i^\(j) without knowing all 
elements of the anomalous dimension matrix ^5J, which was verified by the exact calculations 
in m. 

Using above arguments, we conclude, that at the NNLO level there is only one possible 
candidate for Tuniii)- Namely, it is the most complicated part of the nonsinglet QCD anomalous 
dimension matrix (with the SUSY relation for the QCD color factors Cp = Ca = Nc). Indeed, 
after the diagonalization of the anomalous dimension matrix the eigenvalues 7/(j) and 
in Eq. ((7j) should have this most complicated part as a common contribution because they 
differ each from others only by a shift of the argument (see Eq. ((HJ) and their differences are 
constructed from less complicated terms. The non-diagonal matrix elements of 7a6(j) in Eq. ((Tj) 
contain also only less complicated terms (see, for example, AD exact expressions at LO and 
NLO approximations in Refs. [T] for QCD and ^ for A/" = 4 SYM) and therefore they cannot 
generate the most complicated contributions to 7z(j) and 7z(j). 

Thus, the most complicated part of the nonsinglet NNLO QCD anomalous dimension should 
coincide (up to color factors) with the universal anomalous dimension 'JuniU)- 

NNLO anomalous dimension for J\f = A SYM 

The final three-loop result ^ for the universal anomalous dimension '~funi{j) for A/" = 4 SYM 

■fNote, that in an accordance with Ref. |[7| our normahzation of '^{j) contains the extra factor —1/2 in 
comparison with the standard normahzation (see pQ) and differs by sign in comparison with Vermaseren-Moch- 
Vogt one 
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IS 



where 



7ij) = 7uni{j) = a7Si(j) + a^luIiU) + a^luniU) + « = ^ . (9) 



|7S(j + 2) = -Su (10) 

^7inkj + 2) = (53 + 5_3)-2 5_2,l + 2Si(52 + 5_2), (H) 
SS(j+2) = 2 5_3 52-55-2 5_2 53-3 5_5 + 24 5_2,l,l,l 



+ 6 5-4,1 + 'S'_3,2 + 'S'_2,3 ) ~ 12 5'_3,i_i + 5'_2,1,2 + S -2,2,1 



32 



-i^S2 + 2 Sfj (^3 5_3 + ^3 - 2 5_2,i j - Si (^8 5_4 + 

+ 2 5| + 3^4 - 12 5_3,i - 10 5_2,2 + 16 5_2,i,i) (12) 

and Sa = SaU), Sa,b = Sa,bij), Sa,b,c = Sa,b,c{j) are harmonic sums 

j I j I 

S-U) = H—a^ Sa,b,c,-{j) = E (13) 



m=l ^ ^1 ^ 



s-a{3) = s.^,b,c,...{j) = E^^..-(-), 

m=l '''' m=l 
5_,,,,e,..(j) = i-iy S.a,b,c,..Xj) + S-a,b,c,-{^) {l - i-iy). (14) 

The expression (HH) is defined for all integer values of arguments (see [211 E]) but can be 
easily analytically continued to real and complex j by the method of Refs. |2ZIE1- 

The limit j ^ 1 

The limit j ^ 1 is important for the investigation of the small- a; behavior of parton distri- 
butions (see review and references therein). Especially it became popular recently because 
there are new experimental data at small x produced by the HI and ZEUS collaborations in 
HERA IH21- 

Using asymptotic expressions for harmonic sums at j = 1+u — > 1 (see formulae in Appendix 
at r = —1) we obtain for the = 4 universal anomalous dimension 'Juniij) in E<1- & 

7^(1 + ^) = - + 0{uj'), (15) 
7^(1 + ^) = -32C3 + 0(^'), (16) 
7^(1+^) = 32C3^- 232C4--112OC5 + 256C3C2 + 0(c^') (17) 

00^ (jj ^ ^ 
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in an agreement with the predictions for 7^°\(1 + uj), 'Junii^ + ^) ^-^so for the first term of 
7ini(l + ^) coming from an investigation of BFKL equation at NLO accuracy in jH] ^. 

The limit j — r, r > 

In the case of A/" = 4 SYM the eigenvalue of the BFKL kernel [Sj is analytic in the conformal 
spin \n\, which allows us to continue it to the negative values of |n| [3]. It gives a possibility 
to find the singular contributions to anomalous dimensions of the twist-2 operators not only at 
j = 1 but also at other integer negative points j = 0, —1, —2, ■ ■ -. In the Born approximation 
for the universal anomalous dimension of the supermultiplet of the twist-2 operators one can 
obtain in such way ■juni = 4 a (\&(1) — \E'(j — 1)) 3J, which coincides with the result of the direct 
calculations (see [HUE] and Eq. (|Tn|) ). 

In the NLO approximation the study of the relation between BFKL and DGLAP equation 
at j — r, r > was done in Its extension to the NNLO level is a subject of future 
investigations. Here we present only the AD results in this limit. 

Using formulae presented in Appendix for the harmonic sums calculated in the limit j = 
uo — r ^ — r, r > we can find singularities of 'juniij) © at small u 



,(1) 



/ um \J I 



/ uni \J I 



ijj 



- ^i(r + 1) - (52(r + 1) + (2)^^ + 0{y? 



2^i(r + l) 



1 + (-1)'-)C(2) + 2(-irS2(r + 1) - + OiJ" 



rc(5) 



.(3) 



.(2) 



.(1) 



+ 



+ 



+ 



+ 



UJ 



(18) 

(19) 
(20) 



where 



.(5) 
.(4) 



.(3) 



.(2) 



16(1 + (-1)^), 
-32(l + 2(-l)05i, 

^2 - - 25_2 + 4C2 



f-11 



5^2 - QSl + 4C2 



-4(-l) 





48^_4 - 32^_3Si 







+ 25_3 - - 45'_2,i - 18(2^1 - 5C3 



•i2 
'-2 



(21) 
(22) 

(23) 



(24) 



§ Unfortunately, the results of Refs. contain a misprint. Namely, the coefficient in front of a? obtained 

in the limit j ^ f in Eq. (39) of Ref. '5' should be multiplied by a factor 4. 
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-8^4 - 24(2^-2 + 8(2^2 - 24(3^1 - 3C4 



+( 



85'_4 + 32>S'_3S'i + 325'i5'; 
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2 



+56^4 - 325_3,i - 64Si5_2,i + 32S_2,2 + 32S_25f 
+1652" + 645_2,i,i + 8(2^-2 + 40(2^2 - 32(2^2 - 48(3^1 + 55C4 



(25) 



Here, to shorten notation when presenting resulting expression for singularities of 7„„i(j) we 
dropped argument r + 1 in all harmonic sums. So, the double-logarithmic poles ~ uj~^ and 
~ uj~^ appear in the anomalous dimensions "Juliij) and '~)uli{j)i respectively, only in the case 
of even r values, that is in an agreement with the predictions from the BFKL equation [S]. 

Resummation of 7„„j and the AdS/CFT correspondence 

In the limit j ^ oo the AD results (jl(J|) - (jl2p are simplified significantly. Note, that this limit 
is related to the study of the asymptotics of structure functions and cross-sections at x ^ 1 
corresponding to the quasi-elastic kinematics of the deep-inelastic ep scattering. 

We obtain the following asymptotics for the Af = 4 universal anomalous dimension 'juniij) 
in Eq. Q with 

llU) = -4(lnj+7e)+0(j-'), (26) 
7S.(J) = 8C2(lnj+7e)+12C3 + O(r0' (27) 
7S.(j) = -88C4(lnj+7e) -16C2C3-80C5 + O(j-'). (28) 

Recently there was a great progress in the investigation of the A/" = 4 SYM theory in a 
framework of the AdS/CFT correspondence [2Bj where the strong-coupling limit a^Ac —>■ oo 
is described by a classical supergravity in the anti-de Sitter space AdS^ x S^. In particular, 
a very interesting prediction |29j (see also jHOl) was obtained for the large-j behavior of the 
anomalous dimension for twist-2 operators 

aN 

-f{j) = a{z)lnj, - = 4d (29) 

TT 

in the strong coupling regime (see Ref. jHS] for asymptotic corrections)^ : 

hm a = -z'/' + ^ + (z-'/') . (30) 
z^oo 87r ^ ^ 



^Here we took into account, that in our normalization contains the extra factor —1/2 in comparison 
with that in Ref. 
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On the other hand, all anomalous dimensions 7j(j) and 7j(j) {i = +, 0, — ) coincide at large j 
and our results for ^uni{j) in Eq- © allow one to find three first terms of the small- 2; expansion 
of the coefficient a{z) 

71-2 11 

lim a = -z+ — - — ttV + ... . (31) 
z~*o 12 720 ^ ^ 

For resummation of this series we suggest the following equation for a [3] " 

z = -a + —a^ (32) 
interpolating between its weak-coupling expansion up to NNLO 

d = -z + —z^~ —n^z^ + 0{z^) (33) 
12 72 ^ ^ ^ ' 

and strong-coupling asymptotics 
2^3 ^1/2 , 6 



z 



TT 



+ ^ + {z-'/^) ^ -1.1026 z'/^ + 0.6079 + O [z'^^^) . (34) 



It is remarkable, that the prediction for NNLO based on the above simple equation is valid 
with the accuracy ~ 10%. It means, that this extrapolation seems to be good for all values of 
z. 

Further, for j —>■ 2 let us take into account, that according to the BFKL equation |7j the 
anomalous dimension of twist-2 operators is quantized in the Regge kinematics: 

7 = 1/2 + 21/ + (j - l)/2 = 1 + (j - 2)/2 + tu (35) 

for the principal series of unitary representations of the Mobius group. On the other hand, in 
the diffusion approximation valid near the leading singularity of the t-channel partial wave the 
eigenvalue of the BFKL kernel is 

J-l=uJo-Diy\ (36) 

where and D are the Pomeron intercept and diffusion coefficient, respectively. These quan- 
tities are functions of the coupling constant. We assume, that for the large coupling constant in 
= 4 SUSY the Pomeron intercept approaches the graviton intercept in the AdS^ x 5*5 space 
|35j . which means, that 

Jo = 1 + cuo = 2 - A , (37) 

where A is a small number. Further, due to the energy-momentum conservation (7 = for 
j = 2) the parameters A and D are equal and 7(7) can be expressed near j = 2 only in terms 
of one parameter 

'1 1/A 



7(j) = (j - 2) 



2 1 + Jl + (j - 2)/A, 



(38) 



"Note, that we use the DR-scheme for coupUng constant which removes —1/12 from coefficients of a in 
Eq.(28) of Ref. |1| (see |S1E1I)- 
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The derivative 7'(2) can be calculated from our results in three first orders of the perturba- 
tion theory: 

7'2 = - = z + —z^ z^ + .... 39 

' ^ ^ 2 2A 6 72 540 ^ ' 

Similar to the case of large j for a resummation of this series we used the following equation 
for S = y(2) (see 01) 

— z = -a + -a . (40) 

Its solution at small z is 

a = z^ z-^ + .... (41) 

6 72 432 ^ ' 

One can verify from eqs. (jH^ and (0^, that the prediction for NNLO based on the simple 
equation ()40p is valid with the accuracy ~ 20%. Therefore we can hope, that this method of 
resummation gives us a good estimate also for the behavior of a at large z 



y (2) = 1 - . + 1 - + 1 + (,-1/2) . (42) 

Thus, one obtains for the intercept of the Pomeron in = 4 SUSY from the resummation 
(HUI) at large z the result 

:' = 2-|5,-V.__L,-._o(.-). (43) 

On the other hand, from eqs. (jH^j) and (jHUjl . using also the following relation valid in 
ADS/CFT correspondence for the string energy at j close to 2 [2111201 

E^ = {3 + Vf-A, r = -27, (44) 

we obtain, that the BFKL equation in the diffusion approximation ()35|1 is equivalent to the 
equation for the leading Regge trajectory in the superstring theory 

J = 2 + yt, t = E'/R\ «' = ^A, (45) 

where R is the radius of the anti-de-Sitter space. 

It is naturally to expect that this Regge trajectory remains approximately linear (up to 
corrections to diffusion approximation of the BFKL equation) for all values of t and j. We can 
attempt to use expression (jHH|l also for large z 



7(j)|.^oo = -Vj-2A-V2. (46) 
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This relation is in an agreement with the prediction of A. Polyakov and other authors 

[2Hll2ni**- 

7(j)l(..Hoo = -\e = -^z'/' - + (47) 

providing that 

A = -z~'/\ (48) 

TT 

This number coincides up to 15% with the estimate A = [-\/3/(27r)]z~^/^ obtained in (jlSjl 
from the resummation procedure (j4(jp . We can expect, that expression (j46|) with parameter 
A calculated in (j48j) gives the anomalous dimension of twist-2 operators for z oo and all j 
(neglecting the nonlinearity effects). 

Recently we were informed by Prof. Chung-I Tan, that the above correction A to the 
graviton spin j = 2 coincides in form with that obtained in his unpublished work with R. 
Brower, J. Polchinsky and M. Strassler from the AdS/CFT correspondence. We thank him for 
helpful discussions. 

Conclusion 

Thus, in this paper we constructed the anomalous dimension 'juniij) ^^e J\f = 4 su- 
persymmetric gauge theory in the next-to-next-to-leading approximation and verified its self- 
consistency in the Regge (j 1) and quasi-elastic (j — > oo) regimes. Our result for universal 
anomalous dimension at j = 4 could be used to determine the anomalous dimension of Konishi 
operator jSB] up to 3-loops. It is remarkable, that our results coincide ''"^ with corresponding 
predictions from dilatation operator approach and integrability [27|inH|- The method, devel- 
oped for this construction, can be applied also to less symmetric cases of M = 1, 2 SYM and 
QCD, which are very important for phenomenological applications. For the verification of the 
AdS/CFT correspondence the calculations of the various physical quantities in A/" = 4 SYM 
attract a great interest due to a possibility to develop non-perturbative approaches to QCD. 

We demonstrated above that the expressions interpolating between the week and strong 
regime work remarkably well both in limit j — oo and j — 2. The integrability of the 
evolution equations for the quasi-partonic operators in LLA |1H IT^ is an interesting property 
of A/" = 4 SYM which should be verified on NLO and NNLO level. We hope to discuss these 
problems in our future publications. 
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Appendix A 



For the special functions, 
formulae 



contributing to our anomalous dimensions we have the following 



Sk{j 


-2) = 


S-k{j 


-2) = 


^-2,1,1,1 (i 


-2) = 




-2) = 


^-2,1,1 (j 


-2) = 


^-2,1 (j 


-2) = 




-2) = 


S-2,2,l{j 


-2) = 


^-3,2[J 


- 2) - 


S-2,2{j 


-2) = 




-2) = 


S-2Aj 


-2) = 




-2) = 



--k I ^) ^"nrW 

-ir'l A + C-.(i - (-in + (-1)'^-. 

— l)*^"*"^ — |C4 + 'S'-3,l + S-2,2 — S-4 — S'_2,l,l|, 

+45_4 - 35_3,1 + (1 - (-l)'-)5_2,l,l(00)|, 

1)''+'|-(C2 - ^-2) - Cs - ^C3(l - (-!)'■) + ^-2,1 - 35_3 



33 

C2'S'-2 + 'S'-2,2 + 25'_3,i — 65'_4 + Y5^^ 



_l)-+li|i(5_3 _ 5_ ) + _ 25_3,1 - 5-2,2 + 35_4 1, 

a; la; 4 J 

— 1)'""''^ — |^C4 + S^2,2 — C2S-2 — S^4^, 

-ir+'l[^{s.s + 2(3) + 35_4 - 3(4}, 

-lY+'\-^S.2 - -S-3 - 2-S_4 - S-2,2 + -Cs - 3C4 
+ (l-(-in5_2,2(00)}, 



UJ 



CO cu^ 



I a;"^ a; J 

-iy+'{-S.s + ^ - ^ + C4 + 45_4 - 5-3,1 

la; o;^ a; 

+(i-(-in5_3,i(oo)}, 



(1) 

(2) 

(3) 
(4) 

(5) 

(6) 
(7) 
(8) 
(9) 

(10) 

(11) 
(12) 

(13) 
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where (^k = (1/2'^ ^ — l)Ck and the functions Sa,b,c = Safi^ij + !)• 
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